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Scanning tunneling spectroscopy was performed on single crystals of superconducting 2i7-NbSe2, 
at 300 mK and in a magnetic field, up to 5 T, applied parallel to the a6-plane. This novel field 
geometry allows the quasiparticle density-of-states spectrum to be measured under finite superfluid 
momentum, while avoiding contributions from the vortex-core bound states. At zero field, we 
observed a fully-gapped conductance spectrum with both gap-edge peaks and sub-gap kinks. These 
spectral features show a systematic evolution with the applied field: the kinks close in while the 
peaks move apart in low fields, and the zero-bias conductance has a two-sloped behavior over the 
entire field range, though dipping anomalously at ^ 0.7 T. Our data was analyzed with recent 
theoretical models for quasiparticle tunneling into a current- carrying superconductor, and yielded 
distinct evidence for multiple superconducting gaps coming from various Fermi-surface sheets of 
different topologies, as well as possible implications on the origin of the coexisting charge-density- 
wave order. 

PACS numbers: 71.45.Lr,74.25. Jb,74.55.+v,74.70.Xa 



I. INTRODUCTION 

The layered dichalcogenide 2i^-NbSe2 is a heavily 
studied material, well known for the coexistence of su- 
perconductivity and charge density wave (CDW) order, 
the former occurring below Tc ~ 7.2 K and the latter 
below TcDW ^ 33 K. Recent interest has been in the 
complex band structure of 2i^-NbSe2, particularly in the 
role that it plays in both the superconducting and CDW 
order—. Angle-resolved photoemission (ARPES) and de 
Haas-van Alphen measurements have revealed a Fermi 
surface (FS) with multiple sheets of different topologies 
located at the center and corners of the Brillouin zone^*^. 
Low-temperature ARPES data have shown all four of 
the either 2D or quasi-2D FS sheets to be gapped be- 
low Tc, with various superconducting gap amplitudes and 
anisotropics, while the lone 3D FS sheet at the zone cen- 
ter appears to be ungapped down to 5 K. Temperature- 
dependent ARPES data have also implicated various FS 
regions as possible pockets for CDW nesting, although 
the precise location of these FS pockets and the exact 
mechanism responsible for the CDW order are still un- 
der debate^^. 

Scanning tunnehng spectroscopy (STS) is a power- 
ful microscopic probe for studying both superconducting 
and CDW condensates in complex materials. For 2H- 
NbSe2 , STS has allowed real-space imaging of the CDW 
surface modulations below TcdWi while giving spectro- 
scopic evidence for fully-gapped superconductivity below 
Tc. The application of STS in a c-axis magnetic field has 
also allowed imaging of the superconducting vortex lat- 
tice, by mapping out spatial variations in the tunneling 
density of states (DOS)^. In these earlier STS measure- 
ments, sub-gap structures were seen in the DOS spectra, 
and attributed either to superconducting gap anisotropy 
outside the vortex coi^^ or to pair-potential bound 



states inside the vortex core-^-^. More recent STS mea- 
surements at ultra low temperatures have indicated that 
the sub-gap features outside the vortex cores are due to 
a second superconducting energy gap^^. This interpreta- 
tion suggests that 2i7-NbSe2 is a multiband superconduc- 
tor, with corroborating evidence from heat-capacity and 
thermal-transport measurements^"— as well as muon- 
spin rotation measurements of the vortex-core size^^. 



A timely challenge for the study of 2i^-NbSe2 is to 
reconcile the STS data with the ARPES data. In par- 
ticular, it would be crucial to explain how the multiple 
FS sheets observed by ARPES, showing significant gap 
anisotropy, could give rise to the fully-gapped quasipar- 
ticle DOS spectrum observed by STS. In this paper we 
report new STS measurements taken at 300 mK, where 
we have applied an in-plane magnetic field of up to 5 
T that generates a screening current on the a6-surface. 
This novel STS geometry has the advantage of avoid- 
ing vortex-core states, and allows a measurement of the 
spectral dependence on superfluid momentum without 
compHcations due to radial and angular dependences of 
vortical currents^^^^. We observed a systematic spectral 
evolution with the applied field, with the gap-edge peaks 
shifting to higher energies and sub-gap kinks shifting to 
lower energies in low fields. Furthermore, the zero-bias 
conductance showed a two-slope behavior over the entire 
field range, with a dip anomaly at ~ 0.7 T. Our data was 
analyzed with recent theoretical models for quasiparticle 
tunneling into a current-carrying superconductor to 
yield distinct evidence for multiple superconducting gaps, 
with the pairing involving several FS sheets of various 
topologies. Our results also have possible implications 
on the location of FS pockets responsible for the coexist- 
ing CDW order. 



II. EXPERIMENTAL DETAILS 

The home-built cryomagnetic STS apparatus used in 
this experiment was specially designed for the magnetic 
field to be applied parallel, rather than perpendicular, to 
the sample surface. For tunneling into the c-axis face of 
a superconductor, this field orientation enables a screen- 
ing current to be generated across its a6-surface (see in- 
set of Figure [2]) . The Pt-Ir tips used were conditioned 
using field emission to ensure their robustness, and RF- 
filters were used throughout the wiring to remove recti- 
fication noise. The tip scanning rate was minimized to 
avoid heating due to the piezo motion, and the magnetic 
field was swept slowly to prevent field drifts of the junc- 
tion. The dl/dV conductance spectra were acquired us- 
ing lock-in amplification with the sample biased relative 
to the tip. These experimental refinements enabled us 
to consistently acquire atomically-resolved topographic 
images and also to achieve maximum energy resolution 
down to 300 mK and up to 9 T. The single crystals of 
2i^-NbSe2 used in our experiment were grown by an io- 
dine vapor transport technique^^, and had Tc ~ 7.2 K, 
and upper critical fields ^ 5 T and H^^2 ^ 15 T re- 
spectively for field perpendicular to and parallel to the 
a6-plane. The thick platelet crystals had large a6-plane 
surfaces, which were cleaved just before cooldown. 

III. RESULTS AND DISCUSSION 

Figure [1] shows a dl/dV spectrum taken on a 2H- 
NbSe2 crystal at 300 mK in zero magnetic field. A 
scanning tunneling microscope (STM) topographic image 
showing both the atomic lattice and CDW modulations 
is displayed in the inset. In the spectrum, which was 
normalized by the background dl/dV at 4 mV, a fully- 
developed gap structure is seen, with gap-edge peaks Vi 
at ±1.3 mV and noticeable sub-gap kinks V2 at ^ ±0.7 
mV. Figure [2] shows how the spectrum evolves as a func- 
tion of the applied magnetic field jiioH. As the field is in- 
creased, Vi shift to higher energies, while V2 shift down 
to lower energies, as the zero-bias dl/dV steadily rises 
from zero. Figure [3] more clearly illustrates this spectral 
evolution, by plotting Vi and V2 '^s. field up to 0.2 T, 
when the sub-gap kinks become less discernible. 

To understand the field evolution of our tunneling 
dl/dV spectra, we consider the recent theoretical model 
of Lukic and Nicol^^, which calculates the DOS spec- 
trum for quasiparticle tunneling into a multiband super- 
conductor which is carrying a supercurrent. The super- 
current adds a Doppler term Vp Qs to the quasipar- 
ticle energy dispersion Ek^ where Vi? is the Fermi ve- 
locity and Qs is the superfluid momentum, thus modify- 
ing the quasiparticle tunneling spectrum^^. For a single- 
band 5- wave superconductor, the prior theoretical work 
by Zhang et ali^ has shown that the gap edges in the 
DOS spectrum would first move outward with increasing 
Qs, due to Doppler redistribution of the gap amplitude 
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FIG. 1. Normalized dl/dV spectrum taken on a 2if-NbSe2 
crystal at 300 mK in zero magnetic field, showing a fully- 
developed gap structure with gap-edge peaks at Vi ^ ±1.3 
mV and sub-gap features at V2 ~ ±0.7 mV. The conductance 
at bias between ^ ±0.6 mV is ^ 0, indicating a fully-gapped 
DOS. Inset is an atomically-resolved STM topographic image 
(60 X 60 A^, 1 mV bias, 0.6 nA current, T = 4.2 K), showing 
the CDW modulations, also seen clearly in the corresponding 
Fourier transform (FT) image. 



along the FS. At sufficiently high gr^, finite- momentum 
depairing becomes severe enough that the supercurrent 
density then starts to decrease, causing the gap edges 
to move inward. For a two-band 5- wave superconductor, 
because the two gaps are depaired at different rates by 
the supercurrent, their respective features in the DOS 
spectrum are expected to evolve differently with ^5. In 
our STS data up to 0.2 T, the opposite movements of Vi 
and V2 with field strongly suggest that they are due to 
two different energy gaps. Namely, a small gap which is 
Doppler-breached and thus collapsing at low fields, along 
with a larger gap which persists to higher fields. In fact, 
the early field breaching of the smaller gap is consistent 
with the fast rise of the zero-bias dl/dV with field, as 
seen in Figure [3l 

As further evidence for this two-gap scenario, we also 
analyze how the zero-bias dl/dV evolves over the entire 
field range of our measurement. Figure 2] plots the nor- 
malized zero-bias conductance, Go, measured at 300 mK, 
versus the applied field which was gradually increased to 
5 T. At each field value, measurements were made at 
different positions on the atomically-flat surface and av- 
eraged to yield the mean and the range between maxima 
and minima, indicated by the diamonds and bars respec- 
tively in Figure m A two-slope field dependence is clearly 
seen in this plot, along with an anomalous dip structure 
occurring at ~ 0.7 T, as featured in the inset. Similar 
two-slope field dependences have also been seen for the 
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FIG. 2. Normalized dl/dV spectra taken on a 2if-NbSe2 
crystal at 300 mK, at magnetic fields of up to /jqH = 1.0 
T, with zero-field data for comparison (dashed line). Inset: 
schematic showing the experimental geometry with the STM 
tip tunneling into the c-axis of the sample, while the in-plane 
magnetic field produces a screening current (js) on the ab- 
surface. 



heat-capacity coefficient 7 measured in both 2i^-NbSe2 
and MgB2, and interpreted as distinct evidence of two- 
gap superconductivity in these system^ii^i22. Accord- 
ing to the model of Lukic and Nicol^^, since the zero- 
bias dl/dV is a measure of the quasiparticle DOS at the 
Fermi level, the slope change occurs when the smaller 
of the two gaps begins to close on parts of the FS due 
to finite-momentum depairing. The dip anomaly, on the 
other hand, has also been observed in field-dependent 
thermal-transport measurements in 2i^-NbSe2, though 
its physical origin has not yet been clearly identified^^. 
It should be noted that the dip anomaly appears to co- 
incide with the slope change, thus it is most likely also 
related to the collapse of the smaller gap. It is also worth 
noting that the upper slope in Figured undershoots unity 
when extrapolated to H^^2 ^ 15 T, consistent with there 
being an additional slope jump between 5 T and 15 T, 
as the larger gap also becomes Doppler-breached at the 
higher gr^. 

The two-gap interpretation of our STS data has im- 
portant implications for the scenario of multiband pair- 
ing in 2i^-NbSe2. First, the early collapse of the smaller 
gap with field, well below ~ 15 T, suggests that it 
comes from the 3D pancake-shaped FS sheet around the 
Brillouin zone center F, since an in-plane field cannot 
as effectively induce depairing currents in the other four 
FS sheets, which are barrel-shaped and have relatively 
small dispersions along /c^^. Second, the small size of 
this gap at 300 mK indicates that it should also be weak 
against temperature, thus explaining why ARPES was 
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FIG. 3. Spectral evolution versus magnetic field up to 0.2 T, 
showing the field dependences of the gap-edge peaks Vi and 
sub-gap kinks V2. As the field is increased, Vi increases from 
~ ±1.3 mV to ~ ±1.45 mV, while V2 decreases from ^ ±0.75 
mV to ~ ±0.6 mV and becomes indistinct above ~ 0.2 T. 



not able to resolve any gap on the pancake FS sheet at 5 
K^. Third, although a total of five FS sheets have been 
observed by ARPES in the six-fold symmetric Brillouin 
zone, the two outer FS barrels around F and K are os- 
tensibly 2D and thus cannot robustly couple with the 
tunneling matrix element along kz"^^. Of the remaining 
three FS sheets, one is the 3D pancake already discussed, 
and the other two are quasi-2D barrels with similar super- 
conducting gap maxima, as measured by ARPES"^. This 
variety of FS topologies could plausibly explain why only 
two gap structures are manifested in the c-axis tunneling 
spectrum at 300 mK, even though all five FS sheets are 
evidently involved in the pairing for 2i^-NbSe2. 

Despite this plausible scenario of multiband pairing in 
2i^-NbSe2, there remains a discrepancy between the STS 
and ARPES data over the superconducting gap topol- 
ogy. According to the APRES data of Kiss et a/.^, the 
inner quasi-2D barrel around K has a fairly large gap 
anisotropy, ranging between 1.1 meV and 0.3 meV. Since 
this sheet has finite /c^-dispersion and is thus expected to 
contribute to c-axis tunneling, it is rather puzzling why 
the STS spectrum shows a vanishing DOS up to ~ 0.6 
meV. One possible explanation is that the ARPES data, 
which was taken at 5 K, has enough thermal smearing to 
produce finite low-energy DOS not seen in the STS data, 
which was taken at 300 mK. A more interesting possi- 
bility involves /c-space truncation of the superconducting 
gap anisotropy by the CDW order^^. Namely, pockets on 
the FS which are connected by the CDW wave vectors 
are effectively removed from contributing to the quasi- 
particle DOS spectrum, thus keeping the gap anisotropy 
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FIG. 4. Plot of the normalized zero-bias conductance, Co, 
versus in-plane magnetic field at 300 mK. Data were measured 
at different positions on the surface and averaged to yield the 
mean (diamonds) and the range between maxima and minima 
(bars). A two-slope field dependence is clearly seen in this 
plot, along with an anomalous dip structure occurring at ^ 
0.7 T. 



from fully showing up in STS. For this scenario to explain 
the STS data, the CDW pockets would have to coincide 
with superconducting gap minima on the FS, i.e. located 
at where the K — M line intersects the inner FS bar- 
rels around as in the scheme proposed by Borisenko 
et al.^. Both of these possibilities could explain the ap- 
parent discrepancy between the STS and APRES data. 
In fact, the latter possibility could help to corroborate 
the location of the CDW pockets, which is still under 
serious debate, and thus to identify the mechanism re- 
sponsible for the CDW order in 2i^-NbSe2. To elucidate 



these issues we will carry out a more detailed analysis 
of our STS data, using the multiband tunneling formal- 
ism described above and explicitly accounting for the FS 
topology observed by ARPES. 

IV. CONCLUSION 

In summary, we have performed scanning tunneling 
spectroscopy on single crystals of superconducting 2H- 
NbSe2, at 300 mK and in a magnetic field, up to 5 T, 
applied parallel to the a6-plane. This novel field geom- 
etry effectively avoids vortex-core states, while allowing 
the quasiparticle DOS spectrum to be measured under fi- 
nite superfluid momentum. We observed a fully-gapped 
conductance spectrum showing both gap-edge peaks and 
sub-gap kinks, which evolved systematically with the ap- 
plied field. Our data was analyzed using recent theo- 
retical models for quasiparticle tunneling into a current- 
carrying superconductor, and showed the spectral evolu- 
tion to bear distinct signatures of two superconducting 
gaps. Our results indicate that 2i^-NbSe2 is an osten- 
sibly multiband superconductor, with the pairing pro- 
cess involving several Fermi-surface sheets of different 
topologies. An apparent discrepancy between the STS 
and ARPES data, in relation to the superconducting gap 
topology, remains to be clarified and may have direct im- 
plications on the location of the FS pockets responsible 
for the coexisting CDW order. 
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